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1 Introduction 


For a majority of people, the use of mathematics is limited to counting and 
simple arithmetic, mostly in the handling of money. But since the work of 
Black, Scholes [2] and Merton [7] in 1973, sophisticated mathematics has 
been required in order to conduct transactions in financial instruments such 
as stock options and other kinds of derivative assets. 

A financial asset or security is a claim to some payment. It may take 
the physical form of a piece of paper on which a legal contract specifying 
the claim is written. Assets are frequently traded — bought and sold. We 
are concerned with establishing the monetary value, here and now, of such 
entities. Because if their correct monetary value is unknown, they cannot be 
traded fairly. 

A forward contract is a simple kind of derivative asset. It is an agreement 
to trade a security at a future time T for a price K called the delivery price. 
The party who is going to buy is said to hold the long position in the forward 
contract, and the party who is going to sell holds the short position. 

At the time of writing the forward contract, no money changes hands 
between the two parties, so the initial value of the forward is zero. But at 
some time t between the time of writing the contract and its time of expiry, 
the underlying security value may increase or decrease. So the value of the 
long position will become positive or negative, and the value of the short 
position will, respectively, become negative or positive. 

Suppose that you hold the long position in the contract. At any time 
between the initial writing of the contract and its expiry date, you can go toa 
third party and negotiate a transaction — either selling off your long position 
if its value is now positive, or buying your way out of it if its value is negative. 
Likewise if you are the holder of the short position. Our problem here is to 
determine the correct monetary values involved in such a transaction. 


So we have an underlying security, concerning which the forward contract 
is written, but now both the long and short positions in the forward have 
values (the one is the negative of the other) and are themselves negotiable 
instruments. In other words, the forward contract is a derivative security. 

If the forward contract is written at time 0, and its expiry date is time 
T, let t denote any intermediate time. Let x(t) denote the price at time t of 
the underlying security, and let f(t) denote the corresponding value of the 
forward contract. It is not difficult to express f(t) in terms of other known 
quantities, such as the risk-free interest rate r. It is shown in Hull [3] that 


f(t) = a(t) -— ker (1) 


But it is not so easy to establish the values of other kinds of derivative assets. 
However, Black, Scholes [2], and Merton [7] opened up the subject of option 
pricing. 

The holder of the long position in a forward contract is obliged to purchase 
the underlying asset at time 7’ for the contracted price K. In contrast, the 
holder of the long position in a European call option on the asset has the 
right, but not the obligation to purchase the underlying asset at time 7’ for 
the price K. 

The Black-Scholes formula for the value at time t of a European call 
option is 

a(t)®(d,) — Ke*? 9 8 (dz) (2) 


where 0 < ®(d;) < 1, 7 = 1,2, and ®(d) is a value of a cumulative normal 
probability distribution which depends, through d, on the behaviour of the 
underlying asset price. 

Our purpose now is to examine in more detail the Black-Scholes-Merton 
analysis, and how it may be improved by using the non-absolute, generalised 
Riemann integration of Henstock. 


2 The Black-Scholes Model 


The basic model is Brownian motion. If a process y(t) is a Brownian motion 
(or, in alternative terminology, a Wiener process) it satisfies: 


1. For any t, < ty the increments y(t2) — y(t) are normally distributed; 


2. The sample paths {y(t) : 0 <t < 7} are continuous functions of t. 
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A normally distributed random variable with mean value pz: and variance 0? 
is said to be 
N(u,0°). 


The increments y(t2)—y(t1) of a generalised Brownian motion will have mean 
value and variance each proportional to ty — t;, so the increments are 


N (u(t — t1), 07(t2 — t1)). 


We then say that the process y has a drift rate pp and variance rate 07. The 
Ito calculus [5] can be used to represent y. This approach uses stochastic dif- 
ferential equations (in reality, equations involving stochastic integrals). The 
SDE used to represent a generalised Brownian motion y is 


dy(t) = odw(t) + pdt 


where w(t) is a standard Brownian motion with zero drift and constant unit 
variance rate; so w(t2) — w(t1) is N(0, te — t1) for all ty < te. 

Our model for a security price x(t) is geometric Brownian motion (or 
exponential Wiener process): x(t) is a geometric Brownian motion if Ina(t) = 
y(t) is a Brownian motion. 

If the Brownian motion y(t) has drift rate ys and variance rate o?, then, 
using the Ito calculus [1], we find the following SDE for x(t): 


dx(t) = x(t)(odw(t) + pdt), (3) 


which can be solved using the Ito calculus to give 
1, 
x(t) = exp |ow(t) + (u— 50 ae 


Writing pe — $0” = v, we say that the security price process x(t) has growth 
rate v and volatility o. 

Implicit in this model of the security price process as a geometric Brow- 
nian motion is a probability measure on the process sample space which 
governs the actual changes in time of the security price x(t). But we must 
impose a new and different set of probabilities (in place of the “real-world” 
or actual probabilities which allow for a risk-premium) on the sample space 
in order to accomplish the so-called risk-neutral valuation which is used to 
establish the present value of the derivative price. 


Let P denote the natural or real-world probability measure, as implied 
in the SDE (3) above. Our next step is to impose a different probability 
measure Q on the sample space of paths {2(t):0<t <7}. This Q will be 
the risk-neutral probability measure that we require. 

To see the connection between two probability measures P and Q, we use 
the Radon-Nikodym Theorem. First, we suppose that P and Q are equivalent 
measures; that is, for any event A of the sample space, P(A) is positive if 
and only if Q is positive. 

Let a be the Radon-Nikodym derivative. The Radon-Nikodym Theorem 
tells us that 


This specifies us at time T. Let ¢(t) := EP(B\x(t)). For t > s it can be 


shown that 1 


E® (x(t)|x(s)) = a (C()z(t)|2(s)). 
(The notation ---|xz(s)) indicates that x(t) is determined (or known) for all 
times less than or equal to s, but is unknown or random for all times t greater 
than s.) This establishes the relationship between the “natural” probability 
measure P of the security price process x(t) and any equivalent probability 
measure @ on the spaces determined by selecting values for s and t: 


Vet )ts Seth Oe eta 


We will require to determine such a measure Q which will change a geomet- 
ric Brownian motion with a given, non-zero drift, to a geometric Brownian 
motion with zero drift. 

This is the fundamental step in the valuation of derivatives (and it applies 
equally to the classical theory which we are discussing at present, and to the 
theory based on Henstock integration which we will discuss later). 

In the classical Ito calculus version of derivative pricing the Girsanov 
Theorem is the key to determining Q. (But when we use the Henstock 
integral, a much simpler approach does the trick.) Loosely, the Girsanov 
Theorem says: Jf w(t) is a P-Brownian motion and y(t) is a pre-visible 
process satisfying the boundedness condition 

ae 


EP (os E [sat 
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then there exists a probability measure Q such that Q is equivalent to P, 


ge = exp (- Mo y(t)dw(t) — 4 fo 7(t)*dt) , and 
w(t) = w(t) + fpy(s)ds is a Q- Brownian motion. 


In other words, given the sample space {y} we can impose a probability 
measure on it so that the resulting Brownian motion y(t) has whatever drift 
rate we want. 

As a consequence of this, if x(t) is a geometric Brownian motion with 
SDE 
dx(t) = x(t)(odw(t) + pdt) 


where w(t) is a P-Brownian motion, then we can use the Girsanov Theorem 
with 7(t) = (w—v)/o to obtain a measure Q so that 


dx(t) = x(t)(odtw(t) + vdt) 


where w(t) is a Q-Brownian motion. 

We can now formulate the classical Black-Scholes model. For 0 < t < 
T, let f(t) denote the price process of a derivative security, the underlying 
security price process being x(t). We need an amount of cash which grows 
at the risk-free rate — a cash bond — so we represent this by the bond price 


b(t) = exp|rt]. 


Let o be the volatility of the security and let ys be its growth rate. Define 
the discounted claim process D as follows: 


D(t) = E® (0(T)"'f(D)|ax(t) = B° (exp [-rT] F)|x(), 
where Q has been chosen so that the discounted stock process S, 
S(t) := b(t)~*x(t) = exp (—rt) x(t), 


is a Q-martingale. The construction of D implies D(t) is also a martingale 
under Q. (It should be noted that, where we have used x, in many books S is 
used to represent the stock price process. ) 

In effect, we use the Radon-Nikodym and Girsanov Theorems to select Q 
so that S(t) has zero growth rate. The existence of such a Q is guaranteed, 
in the absence of arbitrage, by the Fundamental Theorem of Asset Pricing 
[4]. 


The next step is to use the Martingale Representation Theorem. As a 
process whose increments are governed by the probability measure Q, S(t) is 
a martingale because of the way Q is selected. D(t) is also a Q-martingale. 
The Martingale Representation Theorem then implies that, if the values of 
S(t) and D(t) are known at time t, there exists a function A(t) so that the 
increments in S(t) and D(t) satisfy the stochastic differential equation 


dD(t) = A(t)dS(t). (4) 


When the value x(t) is known, and hence S(t) and D(t)) are known, the 
function A(t) is deterministic in terms of the SDE just given. But since x 
is stochastic, A(t) is globally stochastic. The term pre-visible is sometimes 
used to describe this property of A. We can think of A as being continuous 
from the left. 

The implication of (4) is that, under the probability measure Q, the value 
of the derivative is growing (“on average”, in terms of expectations relative 
to Q) at the risk-free rate, just as the value of the stock is. 

The basic Black-Scholes result is the following. Suppose r, jz and o are 
constant. Suppose the derivative claim f is determined at time 7. At time 
t, z(s) is taken to be known for 0 < s < t, so x(s) is not a random variable 
for any such s. Since f(t) is determined by x(s) (0 < s < ft), neither f(t) 
nor D(t) are random variables. (Randomness only occurs for times greater 
than t. This is the conditionality implied by the notation ---|z(t).) Then 
the arbitrage price f(t) of such a claim is given by 


b(t) f(t) = ES (oT) f(T) |x(t)) 
f(t) = dU) B? (O(T) f(T) |x(t)) 


This gives us the valuation we need. The derivative value at time t is 


f(t) =e TP VES (F(T) |a(t)). (5) 


iFrom this we can deduce the values of various kinds of derivative instru- 
ments, whose values at maturity are given by appropriate functions f(T). 

This is the outline of the classical Black-Scholes model. The above is, 
for the most part, a summary of an account of the classical continuous-time- 
continuous-values model given in the book [1] of Baxter and Rennie, which 
is a good place to start reading this theory, before attempting the deeper 
mathematical treatments of the subject. 


or 


3. Pricing a European Call Option 


A European call option gives the holder the right, but not the obligation, 
to buy a unit of stock for a pre-determined amount K on a particular date 
in the future, say T. If the stock price x(T7') at time T is less than Kk, the 
option is not exercised and its value is zero. If #(Z’) is greater than K, the 
option is exercised, and a profit of z(T) — K can be realised if the stock is 
immediately sold. 

The value of the claim at expiry time T is therefore 


f(T) = max(x(T) — K,0). (6) 


The Black-Scholes theory then gives the present value of the derivative as 
f(t) = et? YE2 (f(T)|x(t)). Letting t = 0 represent the present, for a 
European call option this reduces to 


f(0) =e E®? (max(x(T) — K,0)), 


where Q is the martingale measure for the discounted stock price process 

S(t). To evaluate f(0) we need an explicit expression for the probability 

distribution of #(T) under Q. If we look at the process for x(t) in terms 

of the standard (drift rate zero, variance rate 1) Brownian motion w(t), we 
have 

d(In x(t)) 

In x(t) 

x(t) 

Thus, if we let z denote a normally distributed random variable with param- 

eters 


odw(t) + (r — $07)dt, 
In x(0) + ow(t) + (r — $07)t, 
(0) exp[ow(t) + (r — $07)t] 


1 
N (—50°7. o°r) ; 
we can write 
a(T) = x(0)e7*"? 


so, letting b = In (5) —rT, 


(O).-= e778 (max((0)e*t"? — K, 0)) 
e* — Ke?) exp (-Se) ds. 


202T 


gaan Ie? (2(0) 


Now let 


s+ s0°T 
v= ; 
oVT 
sO { 
a 1,2 2 
f(0) = — | (al Qe-ir ee se ne) e 2” du 
V 277 J—oo 
where 
= In 20) + (r— $0°)T 
oVT 
Since 
en OV Te 30° T— 50 = eng etovT)? 
we get 
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atovT 1 4 a 1 
f(0) = 2(0) (se I. e7 2” wv) — Ket Ge Re e723” wv) , (7) 


and this is the celebrated Black-Scholes formula for the price of a European 
call option. Compare its form to the form of the expression for the value of a 
forward contract (1) above. If the present value x(0) of the underlying asset 
is substantially larger than the exercise price K in the option, then the option 
is very likely to be exercised, and, to some extent, performs financially like 
a forward contract. In that case a will be large, and the bracketed factors in 
(7) above will approach 1. In another notation, 


f(0) = 2(0)®(a + oVT) — Ke" ®(a) (8) 


where ®(a) is the probability that a standard normal variable v (with mean 
zero and variance 1) is less than or equal to a. 


4 Derivative Pricing and Henstock Integrals 


The classical Black-Scholes-Merton method for pricing European call options 
uses the It6 calculus to model the processes involved. But it is possible to 
model stochastic process using Henstock integrands instead of Ito differentials 
(or stochastic integrals), and to derive the Black-Scholes partial differential 
equation and pricing formulae using elementary methods. This is done in [8]. 
We outline some of the basic ideas here. 


The Black-Scholes model [2] assumes that the price of an economic asset, 
as arandom function of time, is a geometric Brownian motion. This implies 
that if the value x;_; occurs at time ¢;-1, the probability of the outcome 
that, at time t; the process takes a value x; between u,; and v,, is related to 


ae el _ (na; —Inz,;-1) 
207 (t; _ 321) 


he 


where A; is a normalising factor (270°(t; — ft) 

When pricing a derivative asset, such as a European call option whose 
value depends on the movements in the value of an underlying asset, the 
probabilities involved turn out to have the form 


[ swae, (9) 


J 


(10) 

with y being the actual growth rate of the underlying asset values. 
iFrom this, the probability of the outcome that, at times ¢;, the under- 
lying asset price process x takes values x; in the range [u,,v,| for 1 <j <n 
will be given by integrating from u; to vj, 7 = 1, 2,...,n, giving an integral 


of the form 
ite a itl CAG 0| dx ,--+d&p. (11) 


In order to accomplish the risk-neutral valuation described in Section 2 above, 
we construct probability functions which give the underlying asset values x(t) 
a growth rate equal to the risk-free interest rate r. This is done simply by 
replacing the parameter ys by r in (10) above, giving 


2 
alee | = | 1 (72 —Inaj-1 — (r — $0°)(t; =42) | 


G2 t; = tj-1 


tj — =) 


(12) 
which is g;(7), and then forming the density functions 


cs a {Ta o} dx ,--+di&p. (13) 
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These are the risk-neutral probabilities for the asset price process x. That 
the risk-neutral probabilities are, firstly, so easily established, and, secondly, 
specified sufficiently in (13), distinguishes the Henstock integral approach 
sharply from the pricing theory that has developed over the past twenty 
five years or so, as described in Section 2 above. The latter approach re- 
quires that the simple sets of outcomes described above be extended, us- 
ing the Kolmogorov Theorem, to a sigma-algebra of measurable sets in an 
infinite-dimensional sample space whose representative elements are continu- 
ous paths; that the processes involved be represented by appropriate stochas- 
tic differential equations; that a suitable probability measure for the sample 
space be found by means of the Girsanov and Radon-Nikodym Theorems; 
and that the derivative asset valuation be then determined by means of an 
expectation using Lebesgue integration. 

The binomial model of derivative valuation (see [3]), in which only a finite 
number of times t; and a finite number of asset and derivative values are 
considered, is much simpler than the continuous time model in which every 
possible time t¢ is allowed. The difficulty arises because of the complicated 
structures of measurable sets which the continuous-time-continuous-values 
model requires in the classical It6 model described in Section 2 above. 

However, to calculate expectation using Henstock integration, the ma- 
chinery of measurable sets is not required. It is sufficient that the probabil- 
ities be defined for u; < a(t;) < vj, 1 <j <n. And in this situation, it is 
remarkably simple to determine the change of measure needed for risk-neutral 
valuation. The simple sets we have just described are, of course, measurable 
sets in the classical Ito theory of Section 2. But since that theory requires 
us to deal with probability measure on general measurable sets, which are 
harder to visualise than the simple sets we have just described, the classical 
theory uses the very abstract notation and methods of stochastic differential 
equations, and tends to lose sight of the basic transition probabilities of (11) 
which are fundamental to both the Ito and Henstock approaches. 

Let us take, for instance, a European call option, whose claim value de- 
pends on the underlying asset value at time T in a very simple way, as we 
have seen in (6). If we seek to obtain the claim value by trying to compute 
a statistical expectation by integrating the discounted claim value, in n di- 
mensions only, with respect to the probabilities defined by the n-dimensional 
integrals of (11) above, we get a result similar to that which is obtained by 
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the Lebesgue integral-based continuous-time model. To see this, take 
OS to = te Se es tH TS ae Sat, ior 9 = 9; 


and, following the reasoning involved in (5), as discussed in Section 2, esti- 
mate the expectation of 


exp|—rT] max(x(T) — K, 0) 


with respect to the probabilities in (11) above. We may suspect that this 
involves evaluation of an integral of the form 


[of ber den 


(Remember that xo is the present, known value of the underlying asset, and 
is not, therefore, a random variable.) The integrand ¢ is 


e"? max(r, — K,0) x Il gj(r). (14) 


j=l 


When the integration with respect to 7; (1 < j <n) is performed, we get the 
very same expression as that found in Section 3 for the value of a European 
call option. 

Another pointer to a new approach to derivative valuation is obtained 
from the Black-Scholes partial differential equation. The problem of option 
valuation was originally solved (with r and o assumed constant) by Black 
and Scholes [2], not by the the risk-neutral probability method described in 
Sections 2 and 3, but by solving the following partial differential equation, 
where po, to are written as €, T, respectively: 


O 1 OF 
+ as Saath (15) 


OT cae ae 


subject to the boundary condition 
f(T) = max(x(T) — K,0). 


Note that the integrand ¢ of (14) depends on the parameters € = Xo, T = to 
If we perform on @ the various partial differentiations that appear in (15), 


11 


we obtain a partial differential equation which is very similar to the Black- 
Scholes equation. We might then hope to take expectations (using some 
system of integration) involving risk-neutral probabilities, and, switching the 
order of integration and differentiation, obtain the Black-Scholes partial dif- 
ferential equation (15). 

These simplistic observations indicate that there may be a way to avoid 
using the Ito calculus and other advanced mathematical theories in formu- 
lating a continuous time model for pricing derivatives, and this is further 
motivation for examining the problem in terms of Henstock integration in 
the manner of [8]. 
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